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2 1 Abstract 

We calculate open charm and charmonium production in Au + Au reac- 

00 ■ 

tions at y/s = 200 GeV within the hadron-string dynamics (HSD) transport 

t-h ; 

approach employing open charm cross sections from pN and ttN reactions 

m 

that are fitted to results from PYTHIA and scaled in magnitude to the 
available experimental data. Charmonium dissociation with nucleons and 

*o , 

formed mesons to open charm (D + D pairs) is included dynamically. The 



'comover' dissociation cross sections are described by a simple phase-space 
model including a single free parameter, i.e. an interaction strength Mg, 
that is fitted to the J/^f suppression data for Pb + Pb collisions at SPS 
energies. As a novel feature we implement the backward channels for char- 
monium reproduction by DD channels employing detailed balance. From 
our dynamical calculations we find that the charmonium recreation is com- 
parable to the dissociation by 'comoving' mesons. This leads to the final 
result that the total J/^f suppression at y/s = 200 GeV as a function of 
centrality is slightly less than the suppression seen at SPS energies by the 
NA50 Collaboration, where the 'comover' dissociation is substantial and 
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the backward channels play no role. Furthermore, even in case that all di- 
rectly produced J/^> mesons dissociate immediately (or are not formed as a 
mesonic state), a sizeable amount of charmonia is found asymptotically due 
to the D + D — ► J '/ VP + meson channels in central collisions of Au + Au at 
y/s = 200 GeV which, however, is lower than the J/* yield expected from 
binary scaling of pp collisions. 
PACS: 25.75.-q; 13.60.Le; 14.40.Lb; 14.65.Dw 
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I. INTRODUCTION 



The dynamics of ultra-relativistic nucleus-nucleus collisions at SPS and RHIC energies 
are of fundamental interest with respect to the properties of hadronic/partonic systems 
at high energy densities as encountered in the early phase of the 'big bang'. Especially 
the formation of a quark-gluon plasma (QGP) and its transition to interacting hadronic 
matter has motivated a large community for about 20 to 30 years by now [1]. However, 
even after more than a decade of experiments at the Super Proton Synchrotron (SPS) and 
recently at the Relativistic Heavy-Ion Collider (RHIC) the complexity of the dynamics 
has not been unraveled and no conclusive evidence has been obtained for the formation 
of the QGP and/or the properties of the phase transition [2,3] though 'circumstantial 
evidence' has been claimed [4]. 

Apart from the light and strange flavor (u, u, d, d, s, s) quark physics and their hadronic 
bound states in the vacuum (ir,K,<f> etc.) the interest in hadronic states with charm 
flavors (c, c) has been rising additionally in line with the development of new experimental 
facilities. This relates to the charm production cross section in pN, irN, pA and AA 
reactions as well as to their interactions with baryons and mesons which determine their 
properties (spectral functions) in the hadronic medium. 

The charm quark degrees of freedom are of special interest in context with the phase 
transition to the QGP since cc meson states should no longer be formed due to color 
screening [5,6]. However, the suppression of J/\& and mesons in the high density phase 
of nucleus-nucleus collisions at SPS energies [7-11] might also be attributed to inelastic 
comover scattering (cf. [12-19] and Refs. therein) provided that the corresponding J/\&- 
hadron cross sections are in the order of a few mb [20-27]. Theoretical estimates here 
differ by more than an order of magnitude [28] especially with respect to J/^-meson 
scattering such that the question of charmonium suppression is not yet settled. On the 
other hand, at RHIC energies further absorption mechanisms - such as plasma screening 
and gluon scattering - might play a dominant role as suggested in Refs. [29,30] and also 
lead to a substantial reduction of the J/\& formation in central Au + Au collisions. 
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On the other hand, it has been pointed out - within statistical models - that at RHIC 
energies the charmonium formation from open charm + anticharm mesons might become 
essential [31] and even exceed the yield from primary NN collisions [31,32]. However, 
a more schematic model by Ko et al. [33] - including the channels J/\I> + n <-> DD 
suggested that such channels should be still of minor importance at RHIC energies but 
become essential at LHC energies. A similar conclusion has been reached in Ref. [34]. 
One of the prevailing questions thus is, if open charm mesons and charmonia will achieve 
thermal and chemical equilibrium with the light mesons during the nucleus-nucleus re- 
action as suggested/anticipated in Refs. [35-38]. Such issues of equilibration phenomena 
are traditionally examined within nonequilibrium relativistic transport theory [13,39-42]. 

In this work we will calculate open charm and charmonium production at RHIC en- 
ergies within the HSD transport approach [13,16,43] for the overall reaction dynamics 
using parametrizations for the elementary production channels including the charmed 
hadrons D, D, D*, D*, D„ D s , D*, D*, J/tf , V(2S), X2c from NN and ttN collisions. The 
latter parametrizations are fitted to PYTHIA calculations [44] above y/s = 10 GeV and 
extrapolated to the individual thresholds, while the absolute strength of the cross sections 
is fixed by the experimental data as described in Ref. [43]. In the latter work we have 
calculated excitations functions for open charm mesons and charmonia including the J/ty 
suppression by dissociation with baryons and meson ('comovers') using the J/^-meson 
cross sections from Haglin [20]. The centrality dependence for the J/\I> survival proba- 
bility has been presented in Ref. [16] for SPS (y/s = 17.3 GeV) and RHIC energies (y/s 
= 200 GeV), too, for Pb + Pb or Au + Au collisions, respectively. We here extend our 
previous works and include explicitly the backward channels 'charm + anticharm meson 
— > charmonia + meson' employing detailed balance in a more schematic interaction model 
with a single parameter or matrix element |M | 2 , that is fixed by the J/ty suppression 
data from the NA50 collaboration at SPS energies. 

Our work is organized as follows: In Section 2 we will present the results of the 
HSD transport approach for charged hadrons, protons, antiprotons and elliptic flow in 
Au + Au collisions at yfs = 200 GeV in comparison to available data. This presentation 
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is necessary since the open and hidden charm formation and propagation proceeds in a 
dense and hot hadronic environment that should be sufficiently realistic. The elementary 
production cross sections for open charm and charmonia from baryon-baryon (BB) and 
meson-baryon (mB) collisions are presented in Section 3 as well as their interaction cross 
sections with hadrons. A phase-space model will be presented, furthermore, for the char- 
monium + meson dissociation cross sections that allows to implement 'detailed balance' 
for all channels of interest. Section 4 contains the actual calculations for the open and 
hidden charm degrees of freedom for Pb + Pb collisions at y/s = 17.3 GeV and Au + Au 
collisions at y/s = 200 GeV with particular emphasis on the novel aspect, i.e. the char- 
monium reformation by open charm mesons employing 'detailed balance'. A comparison 
to the preliminary data of the PHENIX Collaboration on J/ty suppression in Au + Au 
collisions at y/s = 200 GeV will be presented, too. Section 5 concludes this study with a 
summary and discussion of open problems. 



II. CHARGED HADRONS, BARYONS, ANTIBARYONS AND COLLECTIVE 

FLOW 

Before coming to the actual charmonium and open charm dynamics at RHIC energies 
we have to investigate, if the HSD transport approach based on string, quark, diquark 
(q, q, qq, qq) as well as hadronic degrees of freedom performs reasonably well with respect 
to the abundancy of light hadrons composed of u,d,s quarks 1 . Such a test is essential 
since the dissociation of charmonia on baryon, antibaryons and mesons is directly propor- 
tional to their density in phase space. We recall that in HSD all newly produced hadrons 
have a formation time of tf = 0.8 fm/c in their rest frame and do not interact during 
the 'partonic' propagation. Furthermore, hadronization is inhibited if the energy density 



1 For a more recent survey on hadron rapidity distributions from 2 to 160 A GeV in central 
nucleus-nucleus collisions within the HSD and UrQMD [45] transport approaches we refer the 
reader to Ref. [46]. 
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- in the local rest frame - is above 1 GeV/fm 3 , which roughly corresponds to the energy- 
density for QGP formation in equilibrium at vanishing quark chemical potential n q . Thus 
'hadrons' only exist as quark- ant iquark or quark-diquark pairs at energy densities above 1 
GeV/fm 3 and only can become ordinary hadrons if the system has expanded sufficiently 
We note that this cut on the energy density is the only modification introduced as com- 
pared to the earlier studies in Refs. [16,43,47] and has also been included in the more 
recent systematic analysis in Ref. [46] from SIS to SPS energies. 

In order to demonstrate the applicability of the HSD approach to nucleus-nucleus 
collisions at RHIC energies we show in Fig. 1 the calculated pseudo-rapidity distributions 
of charged hadrons (solid lines) for Au + Au at y/s = 200 GeV for different centrality 
classes in comparison to the experimental data of the PHOBOS Collaboration [48] (full 
points), where the error bars indicate the systematic experimental uncertainty. The open 
squares in the upper left figure correspond to the data from the BRAHMS Collaboration 
for the same centrality class [49]. We find that the HSD calculations show a small dip 
in dN/dr] at midrapidity for all centrality classes, which is not seen in the experimental 
distributions. Furthermore, the pseudo-rapidity distributions are slightly broader than 
the data which also might point towards an improper string fragmentation scheme in the 
LUND model [50] employed in HSD. We expect that this issue can be settled uniquely 
when high statistics data for pp reactions at RHIC energies become available. On the 
other hand, the overall description of the rapidity distributions is reasonable good for our 
present purposes. 

A further question is related to the antibaryon and baryon abundancies at midrapidity 
that show the amount of baryon stopping and antibaryon production [51]. We mention 
that multi-meson fusion channels play a sizeable role in recreating baryon-antibaryon pairs 
[47,52,53] and reducing the number of light mesons accordingly. Thus detailed balance 
on the many-particle level - as only found more recently [47] - leads to an approximate 
chemical equilibrium of antibaryons with mesons whenever the meson density is suffi- 
ciently high as e.g. in nonperipheral Au + Au collisions at RHIC energies. Our numerical 
results for the (p + A)/ (p + A) ratio in 10% central Au + Au collisions at y/s = 200 GeV 
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are displayed in Fig. 2 as a function of rapidity y in comparison to the data from the 
BRAHMS Collaboration [54], that correspond to the measured p/p ratio, however, include 
some still unknown fraction from A and A decays. The comparison in Fig. 2 thus suffers 
from a 5-10% systematic uncertainty. We mention that (within statistics) practically the 
same rapidity distribution for antiprotons is obtained when discarding baryon-antibaryon 
annihilation as well as the backward channels. Thus the calculations for charmonia and 
open charm mesons in Section 4 will be performed in the latter limit. Nevertheless, Fig. 
2 suggests that the antiproton/proton ratio is reasonably described in the HSD approach. 
This also holds for the net proton {p — p) rapidity distribution as seen from Fig. 3 in 
comparison to the preliminary data of the BRAHMS Collaboration [55] for the same event 
class as in Fig. 2 2 . 

In principle, one might argue that a transport approach based on string and hadronic 
degrees of freedom should not be adequate in the initial stage of nucleus-nucleus collisions 
at RHIC energies where a new state of matter, i.e. a quark-gluon plasma (QGP), is 
expected/hoped to be formed. However, the global event characteristics and particle 
abundancies from SIS to RHIC energies are found experimentally to show a rather smooth 
evolution with bombarding energy [56,57] such that no obvious conclusion on the effective 
degrees of freedom in the initial phase can presently be drawn. Moreover, the large 
pressure needed to describe the elliptic flow at RHIC energies is approximately described 
by 'early' hadron formation - as in HSD - and the 'large' hadronic interaction cross 
sections. This is demonstrated in Fig. 4 where we show the calculated elliptic flow t> 2 
for charged hadrons (solid lines) as a function of the pseudorapidity r] (upper part) and 
as a function of the number of 'participating nucleons' N part (lower part) for \r]\ < 1 in 
comparison to the preliminary 'hit-based analysis' data of the PHOBOS Collaboration 
[58]. Note, that the experimental error bars correspond to la statistical errors, only. 

2 The experimental data again include some unknown fraction of A and A decays such that the 
'real' (p — p) rapidity distribution should be slightly lower. 
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Our calculations underestimate the V2(r)) distribution close to midrapidity and also are 
somewhat low in the centrality dependence of the elliptic flow. Whereas the elliptic flow at 
midrapidity is well described by hydrodynamical models, the V2(rj) distribution comes out 
too flat in these calculations [59]. We note, that our HSD results are very similar to those 
of the hadronic rescattering model by Humanic et al. [60,61] and almost quantitatively 
agree with the calculations by Sahu et al. [62] performed within the hadron-string cascade 
model JAM [63]. 

On the other hand, unexpectedly high parton cross sections of ~ 5-6 mb have to 
be assumed in parton cascades [64] in order to reproduce the elliptic flow t>2(j>r) seen 
experimentally. These cross sections are about 1/9 of the baryon-baryon total cross section 
(~ 45 mb) or 1/6 of the meson-baryon cross section (~ 30 mb) such that the effective 
cross section for the constituent quarks and antiquarks is roughly the same in the partonic 
and hadronic phase. In this context it will be important to have precise data on open 
charm and charmonium transverse momentum (p T ) spectra since their slope might give 
information on the pressure generated in a possible partonic phase [65]. This argument 
is expected to hold especially for J/\l> mesons since their elastic rescattering cross section 
with hadrons should be small in the hadronic expansion phase [66]. We note, that in 
central Pb + Pb collisions at SPS energies the spectral slope of J/^/ mesons is found 
experimentally to be substantially smaller (~ 240 MeV [67]) than that of protons (~ 300 
MeV [68]). At RHIC energies the radial flow in central Au + Au collisions is even larger 
leading to a stiffer spectrum with an inverse slope parameter ~ 400 MeV for the strongly 
interacting protons [69]. 

Nevertheless, in addition to nucleus-nucleus collisions from SIS to SPS energies [46] the 
HSD transport approach is found to work reasonably well also at RHIC energies for the 
'soft' hadron abundancies such that the 'hadronic environment' for open charm mesons 
and charmonia should be sufficiently realistic. 
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III. ELEMENTARY CROSS SECTIONS 



In order to examine the dynamics of open charm and charmonium degrees of freedom 
during the formation and expansion phase of the highly excited system one has to know the 
number of initially produced particles with core quarks, i.e. D, D, D* , D*, D S1 D s , D*, D*, 
J/*,*(25), X 2c 

A. Production cross section is pp and ttN collisions 

In Ref. [43] we have fitted the total charmonium cross sections (X = xc, J/^,^') 
from NN collisions as a function of the invariant energy y/s by the function 

<4 N (s) = b x (i-^J (™y e(VS-V^) (i) 

with a — 10, (3 — 1, while y/s~o denotes the threshold in vacuum. The parameters were 
fixed in [43] to describe the J/^f and \&' data at lower energy (y/s < 30 GeV). For our 
present study we use the same parametrization (1) with a slightly modified parameter 
P = 0.775 (instead of f3 — 1) in order to fit the preliminary data point from the PHENIX 
Collaboration [70] at y/s = 200 GeV, which gives a(pp -> J/* + X) = 3.8 ± 0.6(stat.) ± 
1.3(sys.) fib for the total J/\I/ cross section. The parameter bx = 240 Cx nb is proportional 
to the fraction of charmonium states Cx- We choose C Xc = 0.4, Cj/y = 0.46, CV = 0.14 
in line with Ref. [71]. 

For the total charmonium cross sections from nN reactions we adopt the parametriza- 
tion (in line with Ref. [14]): 

<(*)=^(l-^f) 7 (2) 

with 7 = 7.3 and d x = 1360.8 Cx nb, which describes the existing experimental data at 
low y/s reasonably well (cf. Fig. 3 from [43]). 

Apart from the total cross sections, we also need the differential distribution of the 
produced mesons in the transverse momentum px and the rapidity y (or Feynman xf) 
from each individual collision. We recall that x F = p z /p™ ax 2p z /yfs with p z denoting 
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the longitudinal momentum. For the differential distribution in xf from iViV and ixN 
collisions we use the ansatz from the E672/E706 Collaboration [72]: 

dN 

^^:~(1-M) C expHwp T ), (3) 

where b PT = 2.08 GeV -1 and c = a/(l + b/y/s). The parameters a, 6 are choosen as 
ajvTv = 13.5, bpfjy = 24.9 for NN collisions and a n N = 4.11, ft^jy = 10.2 for 7rA^ collisions. 

In Fig. 5 (upper part) we compare the calculated J/4> differential cross section in 
rapidity y cm - multiplied by the branching ratio to dileptons - with the preliminary data 
from the PHENIX Collaboration [70] for pp collisions at y/s = 200 GeV using j3 = 0.775. 
Our elementary J/ty formation is seen to be in sufficient agreement with the preliminary 
data [70] though the rapidity distribution appears slightly broader than the data. 

The number of primary J/ty mesons formed in central Au + Au reactions at y/s = 
200 GeV can be estimated - on the basis of the Glauber model - by multiplying the 
pp production cross section with the number of binary collisions (N bin w 1.2 x 10 3 ) and 
dividing by the inelastic pp cross section (~ 45 mb). This leads to a multiplicity of 
primary J/^'s of ~ 0.1 in very central Au + Au collisions. 

The total and differential cross sections for open charm mesons from pp collisions, 
furthermore, are taken as in Ref. [43]. They also might have to be reduced slightly as the 
charmonia cross sections, however, no experimental constraint is available so far. We thus 
refer to the results of Ref. [43] which give ~16 DD pairs in central Au + Au collisions at 
y/s = 200 GeV, a factor of ~160 relative to the expected primordial J/\I/ multiplicity of 
~ 0.1. Note, that at y/s « 17.3 GeV the primary DD to J/* ratio is about 40 [43]; the 
increase of this ratio by a factor of ~ 4 from y/s = 17.3 GeV to y/s = 200 GeV is within 
the expected range. Our results for the rapidity distribution of open charm mesons from 
pp collisions at y/s = 200 GeV (summing up all D and D mesons) is displayed in the 
lower part of Fig. 5 and shows a rather flat distribution at midrapidity, too. Presently, 
there are no data that could control this open charm rapidity spectrum. 

Apart from primary hard A^A^ collisions the open charm mesons or charmonia may 
also be generated by secondary 'meson'-'baryon' (mB) reactions. Here we include all 
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secondary collisions of mesons with 'baryons' by assuming that the open charm cross 
section (from Section 2 of Ref. [43]) only depends on the invariant energy y/s and not on 
the explicit meson or baryon state. Furthermore, we take into account all interactions 
of 'formed' mesons - after a formation time of t f = 0.8 fm/c (in their rest frame) [74] 
- with baryons or diquarks, respectively. As pointed out in Ref. [43] the production of 
open charm pairs in central Au + Au collisions by mB reactions is expected to be on the 
10% level. 

In order to study the effect of rescattering we tentatively adopt the following dissoci- 
ation cross sections of charmonia with baryons independent on the energy (in line with 
Refs. [16,43]): 

&ccB = 6 mb; oj/^b = 4 mb; a XcB = 5 mb; ovb = 10 mb. (4) 

In (4) the cross section a c5 B stands for a (color dipole) pre-resonance (cc) - baryon cross 
section, since the cc pair produced initially cannot be identified with a particular hadron 
due to the uncertainty relation in energy and time. For the lifetime of the pre-resonance 
cc pair (in it's rest frame) a value of t c5 = 0.3 fm/c is assumed following Ref. [75]. This 
value corresponds to the mass difference of the ty' and J/^> . 

For D,D*,D,D* - meson (ir,r),p, u) scattering we address to the calculations from 
Ref. [22,23] which predict elastic cross sections in the range of 10-20 mb depending on 
the size of the formfactor employed. As a guideline we use a constant cross section of 10 
mb for elastic scattering with mesons and also baryons, although the latter might be even 
higher for very low relative momenta. 

B. Comover dissociation channels 

As already pointed out in the introduction the J/^f formation cross sections by open 
charm mesons or the inverse comover dissociation cross sections are not well known 
and the significance of these channels is discussed controversely in the present litera- 
ture [28,31,32,34,76,77]. Whereas in Refs. [16,43] the energy-dependent J/^-meson cross 
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sections for dissociation to DD have been taken from the calculations of Haglin [20], 
we here introduce a simple 2-body transition model with a single free parameter Mq, 
that allows to implement the backward reactions uniquely by employing detailed balance 
for each individual channel. Since the meson-meson dissociation and backward reactions 
typically occur with low relative momenta ('comovers') it is legitimate to write the cross 
section for the process mi + m 2 — > m 3 + m 4 as 

4 E X E 2 E Z E 4 2 ( M 3 + M 4 \ 6 P f 
^1+2^3+4 ( V s) = 2 \M f \ I — I — , (5) 

where E 1 , and denote the energy and spin of hadron i, respectively. The initial and 
final momenta for fixed invariant energy yfs are given by 
p2 _ (s - (Mi + M 2 ) 2 )(s - (Mi - M 2 f ) p2 _( s -(M 3 + M 4 ) 2 )(s-(M 3 -M 4 ) 2 ) 



4s 1 1 4s 

(6) 

where Mj denotes the mass of hadron i In (5) |M/| 2 stands for the effective matrix 
element squared which for the different 2-body channels is taken of the form 



\M S ? = 


-Ml 


for (7r,p) + J/^ ^ 


+ D (7) 


\M f \ 2 = 


3M 2 


for (tt,p) + J/m ^ 


J D* + J D, J D + J D*, J D* + J D* 


\M f ? = 


\Ml 


for (if, K*) + J/V 


-^D S + D,D S D 


W f \ 2 = 


-Ml 


for (if, K*) + J/tf 





involving a single parameter Mq to be fixed at SPS energies in comparison to the data of 
the NA50 Collaboration [9,10]. The relative factors of 3 in (7) are guided by the sum rule 
studies in [78] which suggest that the cross section is increased whenever a vector meson 
D* or D* appears in the final channel while another factor of 1/3 is introduced for each s 
or s quark involved. The factor ((M 3 + M 4 )/ ^fsf in (5) accounts for the suppression of 
binary channels with increasing y/s and has been fitted to the experimental data for the 
reactions 7r + iV^ p + N, u + N, § + N, K + + A in Ref . [79] . For simplicity we use the 
same matrix elements for the dissociation of Xc and with mesons though there is no 
fundamental reason why these matrix elements should be the same. However, since we 
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here concentrate only on the net J/^f absorption and production and not on the explicit 
charmonium 'chemistry', this approximation should work out reasonably well within the 
range of systematic uncertainties. 

The advantage of the model introduced in (5) is that detailed balance for the binary 
reactions can be employed strictly for each individual channel, i.e. 



and the role of the backward reactions (J/\l/+meson formation by D + D flavor exchange) 
can be explored without introducing any additional parameter once Mq is fixed. The un- 
certainty in the cross sections (5) is of the same order of magnitude as that in Lagrangian 
approaches using e.g. SU(A)fi avor symmetry [22,23] since the formfactors at the vertices 
are essentially unknown [78]. 

As mentioned before, we fit the parameter Mq to the J/\I/ suppression data from the 
NA50 Collaboration for Pb + Pb collisions at 160 A GeV (cf. Section 4.1). For the value 
M5 2 = 0.13 fm/GeV 2 used below we end up with the J/^f dissociation cross sections 



displayed in Fig. 6 with n, p, K and K* mesons. The summation over the final channel c 
in (9) includes all binary channels compatible with charm quark and charge conservation. 
Note, that for the comover absorption scenario essentially the regime 3.8 GeV < y/s < 4.8 
GeV is of relevance (cf. Fig. 7.13 in [13]) where the dissociation cross sections are on the 
level of a few mb. We note, that the explicit channel J/ty + n — > D + D, which has often 
been calculated in the literature [22,23,76,77] is below 0.7 mb in our model. A somewhat 
more essential result is that the J/\& dissociation cross section with p-mesons is in the 
order of 5-7 mb as in the calculations of Haglin [20] used before in Ref. [43], since this 
channel was found to dominate the J/^ dissociation at SPS energies [13]. The explicit 
shape of the cross sections is characterized by a rapid rise in y/s whenever a new channel 
opens up. On the other hand, the channels with vector mesons (p, K*) are 'exothermal' 
and thus divergent at threshold. 



0"3+4-»l+2(v / s) = °"l+2^3+4(\/s) 



(2S 1 + 1)(2S 2 + 1) Pf_ 

(2^3 + 1) (2S-4 + 1) Pf 



(8) 



(V5) = E 



(9) 



c 



13 



The cross sections for the backward channels D + D,D + D*,D* + D, D* + D* -> J/tf 
+ meson as well as the channels involving s or s quarks, i.e. D s + D,D S + D*,D* + 
D,D* + L>+ -> J/* + (K,K*), then are fixed by detailed balance via (8). The actual 
results for these channels - summed up again over all possible binary final states - are 
displayed in Fig. 7 separately for the 'non-strangeness' (upper part) and 'strangeness' 
channels (lower part) showing again divergent cross sections for 'exothermal' channels 
like D + D — > J/ty + n. Such divergent cross sections arise in all 'exothermal' S-wave 
channels implying that D + D or D* + D mesons with low relative momentum have a large 
cross section for c and c quark exchange. In actual transport calculations such divergent 
cross sections impose no problems since the transition rates ~ -P/03+4^1+2 remain finite, 
as it is easily seen when inserting (5) into (8), since the divergent factor P| cancels out. 
Furthermore, in the transport calculations an explicit cut in the total cross sections of 
120 mb is employed, which simulates the screening of large cross sections at finite hadron 
density. 



C. Numerical implementation 

We recall that (as in Refs. [43,73,80,81]) the charm degrees of freedom are treated 
perturbatively and that initial hard processes (such as cc or Drell-Yan production from 
NN collisions) are 'precalculated' to achieve a scaling of the inclusive cross section with 
the number of projectile and target nucleons as Ap x A T when integrating over impact 
parameter. To implement this scaling we separate the production of the hard and soft 
processes: The space-time production vertices of the cc pairs are 'precalculated' in each 
transport run by neglecting the soft processes, i.e. the production of light quarks and 
assosiated mesons, and then reinserted in the dynamical calculation at the proper space- 
time point during the actual calculation that includes all soft processes. As shown in Ref. 
[43] this prescription is very well in line with Glauber calculations for the production of 
hard probes at fixed impact parameter, too. We mention that this 'precalculation' of cc 
production might be modified at RHIC energies due to changes of the gluon structure 
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functions during the heavy-ion reaction or related shadowing phenomena [82]. Such ef- 
fects, however, are expected to be of minor importance at RHIC energies (and below) and 
will be discarded for our present study, that concentrates on the balance between comover 
absorption and J/ty reproduction channels. 

Each open charm meson and charm vector meson is produced in the transport calcu- 
lation with a weight Wi given by the ratio of the actual production cross section divided 
by the inelastic nucleon-nucleon cross section, e.g. 



In the transport simulation we follow the motion of the charmonium pairs or produced 
D, D, D* , D*-mesons within the full background of strings/hadrons by propagating them 
as free particles, i.e. neglecting in-medium potentials, but compute their collisional history 
with baryons and mesons or quarks and diquarks. For reactions with diquarks we use the 
corresponding reaction cross section with baryons multiplied by a factor of 2/3. For 
collisions with quarks (antiquarks) we adopt half of the cross section for collisions with 
mesons. 

Furthermore, in addition to our previous studies [16,43,81] the recreation of charmonia 
by channels such as D* + D — > J/^f + ir etc. is taken into account in each individual run 
according to the cross sections (8) with the weight of the produced charmonium states k 
given by 



where Wi, Wj are the individual weights of the open charm mesons. The open charm 
mesons are not allowed to rescatter within a formation time of 0.3 fm/c (in their rest 
frame) since a finite time is needed to form their wave functions. This formation time is 
not well known and presently can only be estimated. Thus we checked - by performing 
calculations with formation times from 0.3 to 0.6 fm/c - that the physical statements (see 
below) remain robust. As commonly employed in transport simulations, the open charm 
meson pairs, that stem from the same interaction vertex, are not allowed to rescatter with 
each other again unless an intermediate scattering has occurred. 




(10) 



W k = WiWj, 
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IV. NUCLEUS-NUCLEUS COLLISIONS 



A. SPS energies 

We directly step on with the results for the charmonium suppression and start with the 
system Pb + Pb at 160 A GeV to demonstrate that the 'late' comover dissociation model 
(5) is approximately in line with the data of the NA50 Collaboration. The corresponding 
J /iff suppression (in terms of the decay branch relative to the Drell-Yan background 

from 2.9 - 4.5 GeV invariant mass) as a function of the transverse energy Et in Pb + Pb 
collisions at 160 A GeV is shown in Fig. 8. The solid line (HSD'03) stands for the HSD 
result within the comover absorption scenario for the cross sections defined by (5) while 
the various data points reflect the different data releases from the NA50 Collaboration 
[7-10]. Note, that the 2002 data [11] (lower part) no longer indicate the drop at the highest 
E T (for analysis B) in line with the HSD calculations from 1997 [81] and the UrQMD 
results from 1999 [18] (dashed histogram). We mention that the present calculation (solid 
line, HSD'03) agrees with the earlier calculations from Ref. [81] (dotted line, HSD'97) 
very well except for the first E T -bm.. Thus the cross sections presented in Fig. 6 do not 
lead to an overestimation of J /iff suppression at SPS energies. There might be alternative 
explanations for J /iff suppression as discussed in Refs. [14,29,30,73,83] and/or further 
dissociation mechanism not considered here. However, for the purposes of the present 
study it is sufficient to point out that the cross sections displayed in Fig. 6 most likely 
are upper limits. 

In order to provide some information on the relative production and absorption chan- 
nels for charmonia in these reactions we show the calculated J /iff rapidity distributions for 
10% central Pb + Pb collisions at y/s = 17.3 GeV in Fig. 9. The ordering of the different 
lines is as follows: the upper dot-dot-dashed line stands for the rapidity distribution of 
J /ty mesons produced by initial BB collisions while the lowest dot-dashed line reflects 
the rapidity distribution of J /iff mesons from secondary mB collisions that are of minor 
importance at SPS energies. The dashed line corresponds to the J/^'s dissociated by 
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baryons (B); this absorption mechanism is denoted as 'conventional J/^f attenuation' by 
the NA50 Collaboration and also present in p + A reactions. The dotted line ('m abs.') 
gives the rapidity distribution for J/^'s dissociated with mesons ('comover absorption') 
while the full solid line stands for the final J/^ rapidity distribution. 

As mentioned in Section 3, the model (5) allows to calculate the backward channels 
- leading to J/\I/ reformation by open charm + anticharm mesons - without introducing 
any new parameter or assumption. The result for the total J/\l/ comover absorption rate 
(solid histogram) in central Pb + Pb collisions at 160 A GeV is shown in Fig. 10 in 
comparison to the J/^/ reformation rate (dashed histogram) that includes all backward 
channels. Since the rates differ by about 2 orders of magnitude, the backward rate for 
J/^f formation can clearly be neglected at SPS energies even for central Pb + Pb reactions. 
This result is essentially due to the fact that the expected multiplicity of open charm pairs 
is ~ 0.12 in central Pb + Pb collisions at = 17.3 GeV (according to the calculations 
in Ref. [43]). Even in case of 'open charm enhancement' (as suggested in Ref. [84]) by a 
factor ~ 3, where the J/\l/ reformation rate would increase by a factor ~ 9, the backward 
channels still could be neglected. 

Since the 'comover' dissociation cross sections employed should be regarded as upper 
limits, we conclude that no chemical equilibration between mesons, open charm mesons 
and charmonia is achieved dynamically at SPS energies. Note, however, that the trans- 
verse mass Mt spectra for all mesons including open charm and charmonia from central 
Pb + Pb collisions scale according to the HSD calculations (cf. Fig. 18 of Ref. [43]), if 
final state elastic scatterings are omitted. Thus statistical model fits still should work for 
the different hadron abundancies. 

B. RHIC energies 

For central Au + Au collisions at y/s = 200 GeV, however, the multiplicity of open 
charm pairs should be ~ 16, i.e. by about 2 orders of magnitude larger, such that a much 
higher J/\l/ reformation rate (~ N^ E ) is expected at RHIC energies (cf. Ref. [30]). In Fig. 
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11 we display the total J/^f comover absorption rate (solid histogram) in comparison to 
the J/ 1 ® reformation rate (dashed histogram) as a function of time in the center-of-mass 
frame. Contrary to Fig. 10 now the two rates become comparable for t > 4-5 fm/c and 
suggest that at the full RHIC energy of y/s = 200 GeV the J/^f comover dissociation is 
no longer important since the charmonia dissociated in this channel are approximately 
recreated in the backward channels. Accordingly the J/\P dissociation at RHIC should 
be less pronounced than at SPS energies. Moreover, there is even a small excess of 
J/ty formation by D + D reactions in the first 2 fm/c qualitatively in line with AMPT 
calculations by Zhang et al. [85]. 

In order to provide some information on the relative production and absorption chan- 
nels for charmonia in these reactions we show - in analogy to Fig. 9 - the calculated J/^> 
rapidity distributions for 12% central Au + Au collisions at y/s = 200 GeV in the upper 
part of Fig. 12. The ordering of the different lines is as follows: the upper dot-dot-dashed 
line stands for the rapidity distribution of J/^f mesons produced by initial BB collisions 
while the lowest dot-dashed line reflects the rapidity distribution of J/^f mesons from 
secondary mB collisions that are of minor importance also at RHIC energies. The dashed 
line corresponds to the J/^'s dissociated by baryons (B) and corresponds to the 'conven- 
tional J/ 1 ® attenuation'. This distribution is approximately the same as the recreation of 
J/^/'s from D + D annihilation (thin solid line with open circles). The dotted line ('m 
abs.') gives the rapidity distribution for J/^'s dissociated with mesons ('comover absorp- 
tion'); it is slightly lower than the D + D recreation channel. The full solid line stands for 
the final J/^f rapidity distribution which is about a factor of ~ 3 lower than the primary 
production from BB collisions. Since all distributions (within statistics) are practically 
flat for \y cm \ < 2 no strong sensitivity of the J/\l/ survival probability is expected for 
different rapidity cuts in this interval around midrapidity. 

We additionally comment on results of HSD calculations that have been performed 
under the assumption of initial J/ty 'melting' by color screening in a QGP phase as 
advocated in Refs. [5,6]. To this aim we have 'deleted' all charmonia created initially 
from primary BB collisions in the calculation, but evolved the system in time with the 
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same production and absorption cross sections as before. The resulting final J/^f rapidity 
distribution for central Au + Au collisions at ^J~s = 200 GeV is shown in the lower part 
of Fig. 12 by the dashed line in comparison to the final J '/\& rapidity distribution from 
the upper part of the figure (solid line). The comparison demonstrates that even in 
case of complete initial charmonium dissociation a finite amount of J/\l/'s should be seen 
experimentally, which is roughly half of the yield expected from the full calculations and 
essentially due to the D + D production channels. Since the latter cross sections are upper 
estimates, the J/^f yield (dashed line in Fig. 12) also has to be considered as an upper 
limit in this case. 

A note of caution should be added in context with Fig. 12 since the actual rapidity 
distributions might change quantitatively when including a more refined model for the 
matrix elements in (7) especially for the Xc and ty' states. Furthermore, in-medium 
modifications (or self-energy corrections) of the open charm mesons (and charmonia) 
should change the final rapidity distributions to some extent since a lowering of D, D 
masses leads to an increase of J/ty + meson absorption rates and a decrease of the 
backward channel rates [85]. For constant matrix elements as in (7) these modifications 
directly result from an enhanced phase space for absorption and a reduced invariant energy 
for the backward channels. On the other hand, for enhanced D, D masses in the medium 
the J/ty + meson absorption rates will be lowered and the backward channels be enhanced 
accordingly. As argued in Ref. [86] charmonium spectroscopy in p induced reactions on 
nuclei might shed some further light on this presently open issue. Nevertheless, our actual 
results for the J/^f reformation by open charm + anticharm mesons are in qualitative 
and even quantitative agreement with the independent transport studies in Ref. [85] that 
also demonstrate a net reduction of J/ty mesons relative to the extrapolations from pp 
collisions with the number of binary collisions. 

We now turn back again to the HSD results for the full calculations. To quantify the 
final J/\I/ suppression in Au + Au collisions at RHIC we show in Fig. 13 the calculated 
J/ty survival probability Sj/y defined as 
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(12) 



where N^ n and N BB denote the final number of J/ty 



mesons and the number of J/^'s 



produced initially by BB reactions, respectively In Fig. 13 the quantity (12) is displayed 
as a function of the transverse energy Et - in units of the transverse energy at impact 
parameter b = 1 fm - for Au + Au collisions with (solid line) and without inclusion of the 
backward channels (dash-dotted line). In fact, the dash-dotted line is (within statistics) 
identical to the previous calculation in Ref. [16] demonstrating a considerable J/\P 'co- 
mover' suppression for central collisions. When including the reformation channels this 
suppression is substantially reduced and leads to a less effective dissociation of charmonia 
than at SPS energies (middle dashed line). Furthermore, we observe that Sj/y < 1 for 
all centralities and thus no J/\I/ enhancement relative to the primary BB production is 
found from our calculations as claimed in the statistical models of Refs. [31,38,87]. We 
also like to recall that the charmonium 'melting scenario' advocated in Ref. [83] should 
lead to a step-like Et dependence of Sj/q, due to a successive melting of the Xc and J/ty 
and an almost complete disappearance of J/^'s for central collisions. Moreover, as shown 
in Refs. [88,89], statistical models on the partonic or even hadronic level lead to very dif- 
ferent predictions for the J/^f multiplicity as a function of centrality in Au + Au collisions 
at y/s = 200 GeV. Since at RHIC energies the predictions of the 'comover' approach, 
the statistical models and the 'melting scenario' are substantially different, experiment 
should clearly decide about the adequacy of the concepts involved. 

The preliminary data of the PHENIX Collaboration [70] allow for a first glance at 
the situation encountered in Au + Au collisions at yfs = 200 GeV. In order to compare 
with the preliminary data we have performed a rapidity cut \y cm \ < 2 in the calculations. 
In Fig. 14 the J/^f multiplicity per binary collision (times the branching ratio B) is 
shown as a function of the number of participating nucleons A part in comparison to the 
data at midrapidity. Whereas our transport results give a monotonous decrease of the 
J/ty yield (per binary collision) with centrality, the statistical charm coalescence model 
of Gorenstein et al. [89] predicts an increase by about 20% from A part = 100 to 380. Since 
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the statistics (and binning in A part ) is quite limited so far on the experimental side, no 
final conclusion can presently be drawn, however, the data neither suggest a dramatic 
enhancement of J/\& production nor a complete 'melting' of the charmonia in the QGP 
phase. 

V. SUMMARY 

In this work we have performed a first comparison of results from HSD transport calcu- 
lations on meson, baryon, antibaryon production and elliptic flow with the (preliminary) 
data for Au + Au collisions at yfl = 200 GeV from the PHOBOS, BRAHMS and PHENIX 
Collaborations. The HSD transport approach, which is based on quark, diquark, string 
and hadronic degrees of freedom, is found to give a quite reasonable description of the dif- 
ferent observables studied in this work. Only the elliptic flow v 2 is underestimated closer 
to midrapidity - quantitatively in line with the hadron-string cascade calculations in Ref. 
[62] - indicating that there might be 'extra pressure' being generated in the 'prehadronic 
phase'. 

On the other hand, hard probes such as charmonia and open D-meson pairs are 
expected to be sensitive to the initial phase of high energy density where charmonia 
might be 'melting' according to the scenario advocated in Ref. [83], their formation be 
suppressed due to plasma screening [33] or absorbed early by neighboring strings [73]. 
However, charmonia might also be generated in a statistical fashion at the phase boundary 
between the QGP and an interacting hadron gas such that their abundance could be in 
statistical (chemical) equilibrium with the light and strange hadrons [35,87]. The latter 
picture is expected to lead not to a suppression but to an enhancement of J/^f mesons 
at the full RHIC energy if compared to the scaled J/^f multiplicity from pp collisions 
[31]. We recall that the 'hadronic comover' dissociation concept has lead to a ~ 90 % 
J/ty suppression in central Au + Au collisions at y/s [16] due to the high meson densities 
encountered, however, as pointed out in [16], the latter calculations had been performed 
without including the backward D + D — > J/^f + meson channels thus violating 'detailed 
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balance'. 

The focus of this work has been to show the dynamical effects from the backward 
channels for charmonium reproduction by D + D channels employing detailed balance 
on a microscopic level. To this aim we have formulated a simple phase-space model for 
the individual charmonium dissociation channels with a single free parameter Mq (cf. 
Section 3), which we have fixed at SPS energies in comparison to the J /iff suppression 
data of the NA50 Collaboration. In fact, the results for the charmonium suppression are 
practically the same as in the previous HSD transport calculations [16,43,81]. From our 
dynamical calculations we find that the charmonium recreation by the backward channels 
plays no role at SPS energies (cf. Fig. 10), however, becomes substantial in Au + Au 
collisions at y/s = 200 GeV and even is slightly larger than the 'comover' absorption 
channel. This leads to the final result that the total J /iff suppression as a function of 
centrality is less pronounced than at SPS energies, where the backward channels play no 
role. Furthermore, even in case that all directly produced J /iff mesons are not formed as 
a mesonic state (e.g. due to color screening), a sizeable amount of charmonia is found 
asymptotically due to the D + D — > J /iff + meson channels which is almost quantitatively 
in line with the AMPT calculations in Ref . [85] for central Au + Au collisions at \fs = 200 
GeV. Since the cross sections for J /iff + meson absorption employed in this work have to 
be considered as upper limits, the charmonium reformation by D + D — > J /iff + meson 
channels should be lower than the J /iff cross section expected from binary scaling of pp 
reactions. The preliminary data of the PHENIX Collaboration [70] are compatible with 
our full transport calculations (cf. Fig. 14), however, improved statistics and also data 
for light systems such as Ne + Ne and Ag + Ag will be necessary to clarify the issue of 
charmonium suppression experimentally. 
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FIG. 1. The calculated pseudo-rapidity distributions of charged hadrons (solid lines) for 
Au + Au at ^fs = 200 GeV for different centrality classes in comparison to the experimental 
data of the PHOBOS Collaboration [48] (full points), where the error bars indicate the systematic 
experimental uncertainty. The open squares in the upper left figure correspond to the data from 
the BRAHMS Collaboration for the same centrality class [49]. 
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. Au+Au, s =200 GeV, 10% central 




FIG. 2. The (p + A)/(p + A) ratio in 10% central + Au collisions at y/s = 200 GeV as 
a function of rapidity y in comparison to the p/p data from the BRAHMS Collaboration [54]. 
Note, that the experimental data include some unknown fraction of A and A decays such that 
the comparison suffers from a 5-10% systematic uncertainty. 
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FIG. 3. The net proton (p — p) rapidity distribution in central Au + Au collisions at y/s = 
200 GeV in comparison to the preliminary data of the BRAHMS Collaboration [55] for the same 
event class as in Fig. 2. Note, that the experimental data include some unknown fraction of A 
and A decays such that the 'real' {p — p) rapidity distribution should be slightly lower. 
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FIG. 4. The calculated elliptic flow v% for charged hadrons (solid lines) as a function of 
pseudorapidity rj (upper part) and as a function of the number of 'participating nucleons' A part 
for \r]\ < 1 (lower part) for Au + Au collisions at = 200 GeV in comparison to the preliminary 
'hit-based analysis' data of the PHOBOS Collaboration [58]. 
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pp->J/^+X, s =200 GeV 




FIG. 5. The calculated rapidity distribution for J/ty mesons (upper part, multiplied by the 
branching to dileptons) and all open charm mesons (lower part) from pp collisions at ^/s = 200 
GeV in comparison to the preliminary data from the PHENIX Collaboration [70] for J/\I/ + X. 
The D + D pair rapidity distribution is obtained by dividing the result in the lower part by a 
factor of ~ 2. 
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FIG. 6. The J/\I/ dissociation cross sections with tt, p,K and K* mesons as specified in 
Section 3. 
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FIG. 7. The cross sections for the channels D + D,D + D*,D* + D,D* + D* -> J/* + 
meson (upper part) and the channels involving s or s quarks D s + D, D s + Z)*, D* + D, 
D* + Z) + — > J/^f + (X, X*) (lower part) as a function of the invariant energy y/s according to 
the model described in Section 3. 
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FIG. 8. The J/ty suppression (in terms of the fj, + fj,~ decay branch relative to the Drell-Yan 
background from 2.9 - 4.5 GeV invariant mass) as a function of the transverse energy Et in 
Pb + Pb collisions at 160 A GeV. The solid line (HSD'03) stands for the HSD result within the 
'late' comover absorption scenario presented in Section 3 while the dotted line (HSD'97) reflects 
the earlier calculation from Ref. [81]. Upper part: the full dots stand for the NA50 data from 
1995, the full squares for the 1996 data, the open triangles for the 1996 data with minimum bias 
while the open circles represent the 1998 data adopted from Refs. [7-10]. Lower part: the open 
and full symbols indicate the preliminary NA50 data from 2000 (analysis A, B and C) [11]. The 
dashed histogram is the UrQMD result from Ref. [18]. 
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FIG. 9. Calculated J/^f rapidity distributions for central Pb + Pb collisions at ^/s = 17.3 
GeV. The ordering of the different lines is as follows: the upper dot-dot-dashed line stands 
for the rapidity distribution of J '/^f mesons produced by initial BB collisions while the lowest 
dot-dashed line reflects the rapidity distribution of J mesons from mB collisions. The dashed 
line corresponds to the J/^'s dissociated by baryons (B) and the dotted line shows the J/^'s 
dissociated by mesons (m). The full solid line gives the final J /iff rapidity distribution. 
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FIG. 10. The calculated rate of J/ty dissociation reactions with mesons (solid histogram) 
for central Pb + Pb collisions at \fs = 17.3 GeV in comparison the rate of backward reactions 
of open charm pairs to J /\P + meson (dashed histogram) according to the model specified in 
Section 3. 
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FIG. 11. The calculated rate of J/*$> dissociation reactions with mesons (solid histogram) 
for central Au + Au collisions at \fs = 200 GeV in comparison the rate of backward reactions 
of open charm pairs to J/^> + meson (dashed histogram) according to the model specified in 
Section 3. 
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FIG. 12. Calculated J/ty rapidity distributions for 10% central Au + Au collisions at yfs 
= 200 GeV. The ordering of the different lines in the upper part is as follows: the upper 
dot-dot-dashed line stands for the rapidity distribution of J/^S mesons produced by initial BB 
collisions while the lowest dot-dashed line reflects the rapidity distribution of J/*S> mesons from 
mB collisions. The dashed line corresponds to the J/^'s dissociated by baryons (B); this 
distribution is approximately the same as the recreation of J/Ws from D + D annihilation (thin 
solid line with open circles). The dotted line ('m abs.') shows the J/vt's dissociated by mesons 
(m), which is slightly lower than the D + D — > J/^+meson recreation channel. The full solid 
line gives the final J/ty rapidity distribution. Lower part: The solid line is identical to the final 
rapidity distribution from the upper part whereas the dashed line is obtained from HSD 
calculations assuming that all charmonia produced from initial BB collisions are 'melted' in a 
possible QGP phase (see text). 
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FIG. 13. The calculated J/^ survival probability Sj/q, as a function of the transverse energy 
- in units of the transverse energy at impact parameter b = 1 fm - for Au + Au collisions with 
(solid line) and without inclusion of the backward channels (lower dot-dashed line) . The dashed 
line (middle) shows the result from Fig. 8 for the same quantity in Pb + Pb collisions at -y/i = 
17.3 GeV for comparison. 
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FIG. 14. The calculated J/^f multiplicity per binary collision - multiplied by the branching 
to dileptons - as a function of the number of participating nucleons N par t in comparison to the 
preliminary data from the PHENIX Collaboration [70] for Au + Au and pp reactions. 
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